Polaronic Heat Capacity in The Anderson - Hasegawa Model 



Manidipa Mitra, P. A. Sreeram and Sushanta Dattagupta 
S. N. Bose National Centre for Basic Sciences, JD Block, 
Sector III, Salt Lake City, Kolkata-700098, Indi^ 

Abstract 

An exact treatment of the Anderson - Hasegawa two - site model, incorporating the presence of 
superexchange and polarons, is used to compute the heat capacity. The calculated results point to 
the dominance of the lattice contribution, especially in the ferromagnetic regime. This behavior is 
in qualitative agreement with experimental findings. 
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The Anderson-Hasegawa (AH) model is a two-site realization of the basic idea of 

ti 

double - exchange (DE) proposed by Zener !2| almost fifty years ago. In the DE scenario a 



localized spin is visualized to be strongly 'Hund's rule' coupled to an itinerant spin at the 
same site governed by strength Jh-, whereas the itinerant spin can tunnel from site to site 
accompanied by a 'hopping integral' t. Because of large Jh the itinerant spin is polarized 
along the localized spin, and as it hops to a neighboring site, it carries with it the memory of 
its spin polarization, thereby polarizing the neighboring local spin as well. Thus transport is 
correlated with spin ordering of localized moments, leading to concomitant metal-insulator 
transition and magnetic ordering. 

Since its inception the DE concept has undergone several extensions including a superex- 
change process yielding antiferromagnetic coupling between localized moments, as well as 
polaronic modification of hopping. Indeed polaronic contributions are considered to be quite 
important for thermodynamic properties of a doped magnetic system e.g., Lai^xXxMnO^, 
{X = Ba, Ca, Sr etc). Both thermodynamic and transport phenomena in manganites sug- 
est the importance of polaron formation and the consequent localization of charge carriers 
3]. The two - site AH model provides an exactly calculable framework in which some of 
these ideas can be tested, for evaluating measurable properties of manganites in the wider 
context of a lattice. Besides it is important to keep track of the quantum nature of the 
localized spin ^ - for instance, Mn^+ is a spin - | ion - even though in much of the DE 
literature the localized moment is viewed as a classical vector. Such an exact quantum 
treatment of the two-site AH model incorporating the roles of superexchange and polarons, 
and their contributions to phase diagram and heat capacity, are the subject of this Brief 
Report. For manganites superexchange interaction is also influenced by Jahn- Teller (JT) 
coupling j5|], which is however not considered here for the sake of simplicity. 

With the preceding background to the scope and purpose of the present work we start 
from the AH model including superexchange for which the Hamiltonian can be written as 

2 

HdE = -t ^(4^C2r + h.C.) - Si.(fi + JS1.S2. (1) 

T 1 = 1 



2 



Here cj^(cj^) is the creation(annihilation) operator of the itinerant electron at site i having 
spin projection r, Si is the locahzed spin, ffi is the itinerant spin at site i and J is the strength 
of the superexchange interaction between neighboring sites. For our case we consider | | = 
S, i.e. the locahzed spins on all sites are taken to have the same value. 

From the exact eigenvalue of H^e we may take the large Jh limit by expanding upto 
0(1/ Jh) and write an effective Hamiltonian for the two - site one electron case as jy] 

Heff = -^■^§^(4c2 + h.c) + JS,.S2 + A^j(ni + 712). (2) 

Here 5*0 is the magnitude of the total spin (locahzed plus itinerant) given by | S^i + S2 + a |, 
and 

AEj = -'^:^iS' + l), (3) 

where S' = Sq ~ 1/2. The first term in Eq. is the one obtained by Anderson- Hasegawa 
when the localized spin is treated quantum mechanically. The third term, represented by the 
number operators 1^^1(2) for the itinerant electron, modifies the double-exchange mechanism 
in the presence of the superexchange interaction given by the second term in Eq. (0). This 
on-site term, proportional to AEj, we should emphasize, is hitherto not widely considered 
in the literature, and is a direct consequence of the quantum nature of the localized spin. 
The spin index r has been omitted from Eq. (0), for the sake of brevity, as the spin moment 
of the itinerant electron in any case is parallel to the localized moment, in the Jh — ^ 00 
limit. 

We now turn our attention to the polaronic effects. The minimal model which reflects 

lattice carrier interaction on the double-exchange can be introduced by dovetailing the Hol- 

stein mechanism on the Anderson- Hasegawa Hamiltonian. Therefore, in the limit of large 

Hund's rule coupling, we may write a two site, single polaron, Anderson-Hasegawa-Holstein 

Hamiltonian from Eq. as, 

2 2 
H = H^ff + giUJoJ2niibi + bl)+g2UJo ni{b2 + bl)+n2ibi + b\) + uJoJ2^% (4) 

i=l 1=1 

where, gi{g2) denotes the on-site (intersite) electron-phonon coupling strength. Note that 
we have considered a single phonon mode for interatomic vibrations of frequency ujq for 



which bi and are the annihilation and creation operators. The Hamiltonian (jU, without 
the te™ A^.. has been the subjeet of exaet analyt.ea. study f„. S ^ oo and . ^ eases 
and a numerical solution for the ^ = | case 7| . 

We separate out the in-phase mode and the out-of-phase mode by introducing new phonon 
operators a = (61 + b2)/V2 and d = (61 — b2)/\/2 in the Hamiltonian. The in-phase mode 
does not couple to the electronic degrees of freedom whereas the out-of-phase mode does, 
leading to an effective electron - phonon Hamiltonian H^, given by, 

Hd = uoSd + AEj Y^Ui-ti j {c\c2 + h.c.) + g-UJo{ni - n2){d + S) + JS1.S2, (5) 

1=1 ^ ' 

where g- = {gi—g2)/\/2. Following ^ we use a Modified Lang-Firsov (MLF) transformation 
and obtain, 

Hd = e^ifde-^, (6) 

where R = X{ni — 722) (c^^ ~ d), A being a variational parameter related to the displacement 
of the d oscillator. The basis set is given by |±, A^) = ^(4 ± 4) |0)e|A^), where |+) and 
|— ) are the bonding and the antibonding electronic states and |A^) denotes the iVth excited 
oscillator state within the MLF phonon basis. The diagonal part of the Hamiltonian in 
the chosen basis is treated as the unperturbed Hamiltonian {Hq) and the remaining part of 
the Hamiltonian Hi = — Hq, as the perturbation. 

The unperturbed ground state is the |+)|0) state and the unperturbed energy, E^^ = 
— teff + JSi.S2. Where = AEj — ujQ{2g_ — A)A and tg// = t g^qif- exp(— 2A^). 
However, in the exact quantum limit of core spins, for given values of g- and J, E^^^^ can 
have four values corresponding to ferromagnetic (FM), canted 1 (CAl), canted 2 (CA2) 
and antiferromagnetic (AFM) orientation of the two spins for | 5*12 | = | 5*1 + ^2 |= 3,2, 1,0 
respectively. Th parameter A is calculated by minimizing the unperturbed ground state 
energy 3. 

We have evaluated the perturbation correction to the energy upto the sixth order and 
the wave function upto the fifth order. The convergence of the perturbation series is very 
good for t/ti^o < 1. Further, to study the effect of an external magnetic field {h) we include 



a term —gfisiSi + S2)-h to the Hamiltonian in equation (jl)), g being the Lande g factor. 
We assume that the external magnetic field is along the direction of S12 and is expressed in 
units of fieff{= gfiB)=i- 

The ferromagnetic (FM) and antiferromagnetic (AFM) orders are related to S'i2(=| 5"! + 
5*2 I) = 3 and 0, whereas Su = 2 and 1 are referred to as canted 1 (CAl) and canted 2 (CA2) 
states respectively. The Fig. 1 shows the phase diagram for the four possible spin orders for 
our system, in the g^ vs J plane. To study the polaronic character one calculates the static 
correlation functions {niUi)o and (^1^2)0, where Ui and U2 are the lattice deformations at 
sites 1 and 2 respectively, produced by an electron at site 1 j8j]. The locations of the large 
polaron region (A) and the small polaron region (B) are indicated in the g^ vs J phase 
digram (Fig. 1). Different ground states, required for our calculation below, can be located 
from the phase diagram in Fig. 1. As our phase diagrams are very similar to the ones 
recently presented by Capone and Ciuchi ^ we henceforth focus only on our new results 
for the heat capacity. 

As mentioned earlier our main emphasis in this Report is on heat capacity based on Eq. 
(jSl). Recently, there have been many specific heat measurements of the colossal magnetore- 
sistance (CMR ) ra anganites at low temperatures with and without an external magnetic 
field y, Q, 11, 3|- According to experiments the specific heat Cy has contributions from 

Ifll 12 1, of many 



conduction electrons, lattice and spin waves. The low temperature data 
CMR materials, show a temperature dependence of the form Cy = 'jT + PT^ + ST^^"^, here 
7, P and 6 are constants. The term 7 arises from charge carriers and it is proportional to 
the density of states at the Fermi level and (3T^ is associated with the lattice contribution, 
P being related to the Debye temperature. The term 5T^^'^ gives the spin wave contribution, 
where the coefficient 6 governs the spin wave stiffness. Okuda et al [Ql have estimated the 
electronic specific heat for Lai_xSrj.Mn03 in the ferromagnetic regime and concluded that 
the carrier mass-renormalization near the metal-insulator transition at x = 0.16 is minimal. 
They have also observed a decrease in the low temperature Cy in the presence of a magnetic 
field. Motivated by these observations, we have carried out a calculation of the specific heat. 
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based on the partition function of the system which, from a cumulant expansion upto the 
2nd order, is given by isl - 



Z{(3) = Zo{f3)exp{ f d^' f d(3"{H,{(3')H,{n)), (7) 
Jo Jo 

where Zo{(3) = Tr(e-^^o) ; Hi{(3) = e'^^"Hie-'^^'\ and (3 = The expression () denotes 

the usual canonical averaging. The specific heat is then calculated (in arbitrary units) from 
the well known relation: 

In the low temperature regime only the zero-and one-phonon states contribute. 

As the specific heat has a bearing on fundamental properties of CMR materials it is 
important to address whether the core spins should be treated classically {S oo) or 
quantum mechanically {S = |) for its theoretical estimation. The difference in the quantum 
and classical cases for specific heat, as far as the core spins are concerned, is exemplified in 
Fig. 2(a) and Fig. 2(b) for FM and AFM cases respectively. Temperature is expressed as 
T = kBTujQ. The quantum case only allows for discrete values of the relative angle between 
two core spins, while in the classical case the angle varies continuously. The quantum results 
evidently yield the correct zero temperature limit. 

In the two-site single polaron model we do not have any scope to vary the carrier concen- 
tration and probe different magnetic states. But we can identify the FM and AFM states 
in vs J phase diagram, in which the FM state is stable for lower J values and the AFM 
ground state is obtained for larger J values. In Fig. 3 we show the CyjT vs curves 
for the FM and AFM cases. It is evident that in the FM case behavior of specific heat 
is more pronounced which is in qualitative agreement with the results of jlfll . As the 
stiffness coefficient 5 is proportional to J [14], the spin wave contribution is not significant 
in FM limits. However in the AFM limit (J = 0.2) the variation of CyjT deviates from the 

law and the spin wave contribution is non negligible. These findings are in qualitative 
agreement with the measurements of Smolyaninova et al Moreover the magnitude of 



,0. 



Cy jT for the AFM limit is higher than that in the FM hmit, as in Smolyaninova et al| 
The suppression of Cy in the FM regime can be intuitively ascribed to the absence of spin 
wave fluctuations, as mentioned earlier. However, while calculating heat capacity in 5* — oo 
limit it can be shown that spin wave contribution (i.e. Cy oc T^/^ ) is dominant in FM case. 
This is due to averaging over all possible relative orientations of core spins. 

We show in Fig. 4 the variation of the specific heat in the low temperature region in 
the FM state with zero and one phonon states. With application of an external magnetic 
field /I, Cy takes lower values than for /I = which is expected, as the average energy 
decreases with apphcation of h in the FM state. For CA1(| S^i2 |= 2), CA2(| S\2 |= 1) 
and AFM (| 5*12 |= 0) states the external magnetic field will tend to align the core spins 
to ferromagnetic order(| Sn |= 3). For CAl, CA2 and AFM states at low field and low 
temperatures it can be shown from the present calculation that Cy increases from the h = 
limit as long as h does not shift | 5*12 | to higher values. For larger h, as the ground state 
changes from lower | S12 \ to higher ones, Cy decreases in the low temperature region. For 
CMR materials, there are some reports on measurements of field dependence of Cy in the 



FM state 



and also in the half doped case It was found that for low doping regions. 



Cy decreases with an increasing magnetic field 9| . However the half doped material showed 
Cy /T as independent of applied field ^|. The present calculation of the external magnetic 
field dependence of Cy qualitatively agrees with these experimental findings in the FM limit. 

In conclusion, the present calculation of Cy using an exactly solvable model reveals some 
of the important features of the double exchange polaronic system. The discreteness associ- 
ated with the effective hopping as a result of the quantum nature of the local spin was shown 
to have a significant consequence for thermodynamic properties. As analytic calculations of 
the heat capacity for CMR material to fit experimental results are not starightforward, be- 
cause of the involvement of several parameters, the present calculation for a simplified model 
indeed serves an important role in indicating general trends. Further, a comparison of the 
computed Cy with measured values underscores the importance of the quantum nature of 
the local spin, a fact often ignored in the current CMR literature. 
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Figure Captions : 



FIG. 1. The g_ vs J phase diagram {h = 0) for | Si \ = \ S2 |= |, t = 1 and cuq = 1. (A) 



and (B) denote large polaron and small polaron region respectively. 

FIG. 2. Variations of Cy (in arbitrary units) with temperature T(= fc^TciJo) for h = 0, 
t = 1, cuo = 1, in classical (sohd line) and quantum (dashed line) formulations of the core 
spins for (a) FM ground state(5'_ = 0.2, J = 0.02) and (b) AFM ground state ((?_ = 0.6, 
J = 0.2). 

FIG. 3. Cv/f vs f 2 for /i = 0, t = 1, Wo = 1 in (a) FM ground state and (b) AFM 
ground state for facilitating comparison with experiments. Here Cy is in arbitrary unit and 
T = kftTujn. The results of Fig. 3(a) are in qualitative agreement with Fig. 2 of Hamilton 



FIG. 4. Variations of Cy (in arbitrary units) for g_ = 0.6, J = 0.01 and t = 1, ujq = 1, 
for different values of the magnetic field h = 0, 0.01, 0.05, which exemplify the magnetic 
field dependence of Cy- 



et al 
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